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method  between  a depth  of  2U  and  35  ft. 

Assuming  a peak  acceleration  of  0.2  g,  the  Seed-Idriss  Simplified  Procedure 
predicts  liquefaction  below  a depth  of  25  ft.  The  empirical  method  predicts 
liquefactions  between  a depth  of  25  and  60  ft  and  between  a depth  of  85  and 
115  ft. 

The  soils  below  the  reactor  at  the  LaCross  site  are  predicted  to  strain 
badly,  if  an  SSE  which  produces  0.12  g at  the  ground  surface  occurs;  and  are 
predicted  to  experience  excessive  strains  and  liquefaction  if  the  SSE  with  n 
peak  acceleration  of  0.2  g occurs.  Because  of  the  limitations  in  the  current 
state  of  knowledge  concerning  liquefaction  and  because  of  the  limited  data  for 
use  in  this  analysis,  WES  cannot  conclude  that  the  reactor  vessel  foundation 
is  safe  if  the  0.12  g SSE  occurs  and  concludes  that  the  reactor  vessel  founda- 
tion is  unsafe  if  the  0.2  g SSE  occurs. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con- 
verted to  metric  (Si)  units  as  follows: 


Multiply  To  Obtain 


inches 

25.  U 

millimetres 

feet 

0.301*8 

metres 

pounds  (mass) 

0.1*535921* 

kilograms 

pounds  (mass)  per 
cubic  foot 

16.0181*6 

kilograms  per  cubic 
metre 

pounds  (force)  per 
square  inch 

6891*.  757 

pascals 

pounds  (force)  per 
square  foot 

1*7.88026 

pascals 
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LIQUEFACTION  ANALYSIS  FOR  LACROSS  NUCLEAR  FOWEH  STATION 


Introduction 

Background 

1.  The  United  States  Nuclear  Regulatory  Commission  (NRC)  requested 
that  the  Waterways  Experiment  Station  (WES ) review  certain  foundation 
conditions  at  the  LaCross  Nuclear  Station,  an  operating  nuclear  power 
plant.  This  plant,  which  is  located  near  LaCross,  Wisconsin,  is  among 
the  oldest  in  the  country  and  was  put  into  operation  before  the  present 
Site  Analysis  Report  review  system  came  into  effect.  This  report  documents 
WES'  review  of  the  LaCross  Nuclear  Station.  Specifically,  the  question 
examined  was  the  earthquake  safety  of  the  pile  foundation  which  supports 
the  containment  vessel.  The  piles  are  driven  through  low  to  medium  rela- 
tive density  sands  and  terminate  in  a dense  sand  layer  approximately  28  ft* 
above  bedrock. 

Scope  of  work 

2.  The  investigation  of  the  foundation  at  the  LaCross  Nuclear  Station 
included  the  following: 

a.  Review  of  Chapter  3,  Soil  Engineering  Properties  contained  in 
the  Application  for  Operating  License  for  the  LaCross  Boiling  Water  Reactor 
by  Pairyland  Power  Cooperative  including  portions  of  Appendix  A,  entitled 
"Field  Exploration  and  Laboratory  Tests, and  associated  design  drawings. 

b.  The  performance  of  a liquefaction  analysis  using  the  Seed- 

2 . . . . , 

Idriss  Simplified  Procedure  assuming  an  earthquake  with  a peak  accelera- 
tion of  0.12  and  0.2  g,  respectively. 

* A table  for  converting  U.  S.  customary  to  metric  (SI)  units  is  given 
on  page  3. 


U 


c.  The  performance  of  a liquefaction  analysis  using  Seed's 
empirical  method  assuming  hoth  the  0.12  and  0.2  g earthquakes  and  compari- 
son with  a "rule  of  thumb"  based  on  the  Japanese  experience  at  Niigata 
in  1961*. 

3.  The  objective  of  this  study  was  to  evaluate  to  the  degree  possible 
with  the  data  available  from  prior  field  and  laboratory  studies  by  others, 
the  seismic  stability  of  the  pile  foundation  which  supports  the  contain- 
ment vessel  at  the  LaCross  Nuclear  Station. 

Review  of  Previous  Work 

U.  As  stated  earlier,  portions  of  Reference  1 were  reviewed  to 
determine  the  soil  profile  under  the  containment  vessel,  including  soil 
properties.  Logs  of  Borings  B-3  and  B-l*  drilled  by  Raymond  International 
in  July  1962  and  Borings  DM-1  and  DM-3,  drilled  under  Dames  & Moore's 
supervision  in  1973,  were  reviewed.  Figure  1 shows  an  idealized  soil 
profile  in  the  vicinity  of  the  reactor  building.  The  ground  surface  is 
at  elevation  (el)  636  ft  mean  sea  level  (msl).  The  groundwater  table 
was  assumed  at  a depth  of  13  ft.  Top  of  bedrock  is  located  at  a depth 
of  about  133  ft. 

5.  During  construction,  the  soil  was  excavated  to  el  615  ft  msl,  and 
piles  were  driven  below  the  reactor  containment  vessel.  These  piles 
terminated  at  el  535  ft.  The  piles  were  50-ton  cast  in-place  concrete 
piles  with  a tip  diameter  of  8 in. , a butt  diameter  of  12  in.  and  an 
outer  shell  of  7-guage  steel  monotube,  and  were  driven  approximately 
3-1/2  ft  on  centers.  No  mention  was  made  of  any  internal  reinforcing 
steel  in  the  piles  in  the  plans  and  reports  provided  to  WES.  A total  of 
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approximately  230  piles  were  driven.  The  data  provided  indicate  that 
the  hammer  used  was  a McKierman-Terry  C-5  double-acting  hammer  with  a 
rated  striking  energy  of  16,000  ft-lb  per  blow.  The  piles  were  driven 
to  at  least  a resistance  of  6 blovs  per  in.  for  the  final  2 to  3 in. 

The  number  of  blows  in  the  last  foot  actually  ranged  from  75  to  330. 

6.  The  soil  below  the  reactor  building  (el  615)  consisted  of  a 
fine  to  medium  sand  with  occasional  zones  of  clayey  silt,  coarse  sand, 
and  fine  gravel,  down  to  an  elevation  of  approximately  535  ft.  At 

el  535  ft,  a 10- ft -thick  fine  to  medium  sand,  with  fine  to  medium  gravels 
is  encountered.  Below  this  gravelly  sandy  layer,  is  an  l8-ft-thick  layer 
of  sand  which  immediately  overlies  the  bedrock. 

7.  Also  show  on  Figure  1 are  the  average  blow  counts,  water 
content,  dry  and  wet  density,  and  shear-wave  velocities  for  the  six 
layers  in  the  idealized  soil  profile.  The  reader  is  cautioned  that  the 
blow  counts  may  or  may  not  be  Standard  Penetration  Test  (SPT)  N values. 

It  was  not  explicitly  stated  in  the  available  source  of  information  how 
the  penetration  tests  were  conducted  nor  were  they  called  "Standard 
Penetration  Test"  results.  The  values  shown  on  Figure  1 are  considered 
approximate  average  values  for  the  layer.  Both  the  blows  per  foot  and 
tne  dry  density  values  were  obtained  from  an  evaluation  of  the  boring 
log  data  in  Reference  1.  The  wet  densities  and  shear-wave  velocities 
which  are  shown  on  Figure  1 are  estimates  based  on  WES'  experience  and 
data  presented  in  Reference  1. 

8.  Figure  2 is  a plot  of  blovs  per  foot  versus  depth.  The  data 
obtained  in  Boring  B-3  are  believed  to  have  been  obtained  prior  to  pile 
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driving.  The  data  obtained  in  Borings  DM-1  and  DM-2  are  believed  to  have 
been  obtained  after  pile  driving.  The  reader  is  cautioned  that  a 1-to-l 
comparison  is  impossible  because  the  data  obtained  in  Boring  B— 3 were 
obtained  using  a 2-in.  split-spoon  sampler  while  the  data  obtained  in 
Borings  DM-1  and  DM-3  were  obtained  using  the  Dames  & Moore  sampler  which 
is  3-1 /l»  in.  in  diameter.  Figure  3 is  the  dry  density  information 
obtained  from  samples  obtained  from  Borings  DM-1  and  DM-3. 

9.  Figure  U is  a plot  of  overburden  pressure  versus  depth  for  th* 
site.  As  stated  previously,  the  water  table  was  assumed  at  a depth  of 
13  ft.  Below  this  depth  both  total  and  effective  overburden  pressures 
are  shown. 

10.  The  blow  count  values  were  assumed  to  be  Standard  Penetration 


Teat  N values  and  were  used  to  compute  relative  density  from  the  following 
equation: 1 


11.7  ♦ 0.76 


222(N)  + 1600  - 53(ao)  - 50(C^)‘ 


1/2 


(1) 


where 


D » Relative  density 
r 

N * Standard  Penetration  Test  N values 

a * Effective  overburden  pressure  in  psi 
o 

Cu  - The  coefficient  of  uniformity 

Using  equation  1,  the  relative  density  of  the  top  105  ft  is  predicted  to 
be  between  50  and  60  percent. 

11.  Review  of  the  available  data1  indicate  that  the  material  has 
the  minimum  density  of  about  100  pcf  and  a maximum  density  of  about  120  pcf. 
These  tests  were  run  on  bulk  samples  obtained  by  combining  representative 
materials  encountered  at  the  site.  WEB*  experience  indicates  that  when 
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materials  are  combined  the  maximum  and  minimum  density  generally  increase. 
This  increase  can  be  as  much  as  8 pcf.  Consequently,  WES  believes  that  a 
minimum  density  of  92  pcf  and  a maximum  density  of  112  pcf  may  be  more 
realistic  for  the  in  situ  material.  The  relative  densities  predicted  by 
equation  1 appear  to  be  more  nearly  the  same  as  those  which  would  be 
obtained  using  the  WES'  maximum  and  minimum  density  estimates  and  the 
in  situ  dry  unit  weights  given  in  Figure  3.  An  analysis  of  pile  geometry 
records  supplied  by  NRC  suggest  an  average  density  increase  of  approxi- 
mately 1 pcf  due  to  pile  driving.  This  is  based  on  the  reduction  of 
void  ratio  which  would  occur  assuming  the  soil  displaced  by  the  pile  went 
entirely  into  taking  up  the  voids  of  the  adjacent  soil.  This  assumes  no 
soil  heave  and  does  not  account  for  any  densification  due  to  vibrations 
during  driving.  An  increase  of  1 pcf  is  not  significant  and  is  believed 
not  to  contribute  substantially  to  the  stability  of  the  soils.  Records 
of  ground  surface  movement  during  pile  driving  were  sought  but  no  such 
information  was  provided.  While  it  is  possible  that  more  densification 
may  have  occurred,  there  were  no  data  made  available  to  WES  which  would 
support  this  hypothesis. 

12.  Dames  & Moore  determined  the  liquefaction  potential  of  the 
subsurface  soil  by  performing  11  stress-controlled  dynamic  triaxial 
compression  tests  on  representative  samples  of  the  material  considered 
to  be  potentially  susceptible  to  liquefaction  during  the  SSE.  Eight  of 
the  samples  were  reconstituted  to  approximately  what  Dames  & Moore  believed 
to  be  the  in  situ  density.  In  addition,  three  reconstituted  samples  were 
tested  at  slightly  greater  densities  than  the  average  in  situ  value  (as 
determined  by  Dames  & Moore)  in  order  to  examine  the  influence  of  density 
variation  on  liquefaction  potential.  These  tests  were  run  at  confining 
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pressures  of  1000  and  2000  psf  and  the  results  of  these  triaxial  tests 
are  shovn  on  Figures  5 and  6.  The  data  points  identified  hy  10b,  111, 
and  112  (the  dry  density  in  pcf)  on  Figure  5 were  for  specimens  consoli- 
dated to  an  effective  confining  pressure  of  1000  psf.  The  data  identified 
by  106  and  105  (dry  density  in  pcf)  on  Figure  6 were  obtained  by  consoli- 
dating the  specimens  to  an  effective  confining  pressure  of  2000  psf. 


lamic  Strength  of  the  Soil 


13.  In  order  to  evaluate  the  liquefaction  potential  of  the  soil 
in  question,  the  laboratory  cyclic  triaxial  test  results  obtained  by 
Dames  & Moore  were  used.  Figures  5 and  6 are  reproductions  of  Dames  & 
Moore's  test  results.  The  as-tested  dry  densities  of  the  remolded  speci- 
mens made  from  material  taken  from  Boring  5 at  a depth  of  about  8-1/2  ft 
are  shown  next  to  the  data  points  for  these  tests  on  Figure  5*  A curve 
has  been  drawn  through  the  data  points  for  a dry  density  of  about  111  pcf. 
As  stated  previously,  WES  believes  that  this  material  is  at  an  in  situ 
dry  density  of  about  102  pcf.  Consequently,  a curve  more  or  less  parallel 
to  the  111-lb  curve  was  drawn  through  a data  point  at  10k  pcf.  This 
curve  (marked  y^  = 102  pcf)  was  used  to  evaluate  the  soil  strength.  On 
Figure  6,  the  remolded  soil  specimens  taken  from  Boring  3 at  a depth  of 
35.5  ft  have  been  used.  Adjacent  to  each  data  point,  the  as-tested 
density  is  listed.  A curve  is  drawn  through  these  data  for  a density  of 
about  106  pcf.  A curve  more  or  less  parallel  to  this  curve  has  been 

drawn  and  labeled  y.  = 100  pcf,  because  WES  believes  the  in  situ  density 

d 

of  the  material  at  a depth  of  about  38.5  is  100  pcf.  This  curve  was 


also  used  to  evaluate  the  dynamic  soil  strength. 
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lU . As  will  be  mentioned  later  on  In  this  report,  the  number  of 

equivalent  cycles  of  load  for  SSE  was  chosen  to  be  10  (see  discussion 

of  design  earthquake  on  page  8 ) . Figures  5 and  6 were  entered  at 

10  cycles  and  the  stress  ratio  required  to  cause  10  percent  double-amplitude 

strain  was  determined.  This  stress  ratio  was  multiplied  by  1000  and 

2000  psf  as  appropriate  to  determine  the  superimposed  dynamic  shear 

strength.  These  values  are  plotted  on  Figure  7«  A line  was  drawn  through 

these  data  points  and  labeled  isotropic  laboratory  data.  Isotropically 

consolidated  cyclic  triaxial  test  data  must  be  corrected  by  a correction 

factor,  Cy  , to  represent  field  conditions.  This  correction  factor  is 

based  on  the  comparison  of  cyclic,  triaxial  test  results  to  cyclic 

simple  shear  and  SHAKE  table  test  results  (References  2,  1* , and  5)«  For  this 

investigation,  a C of  0.57  was  used.  Also  shown  on  Figure  7 is  a 
Y 

curve  labeled  field  conditions.  The  ordinate  of  this  curve  is  0.57 
times  the  ordinate  for  the  curve  labeled  isotropic  laboratory  data.  The 
field  condition  curve  was  used  to  evaluate  the  dynamic  shear  strength  of 
the  soil  during  this  investigation. 

Seed-Idriss  Simplified  Procedure 

15.  In  order  to  evaluate  the  liquefaction  potential  of  a site, 
the  cyclic  stresses  generated  by  the  earthquake  must  be  determined. 

Reference  2 suggests  that  the  average  shear  stress,  TRve  , generated 
by  an  earthquake  can  be  determined  by  the  formula: 


T 


ave 


0.65 


A . x 


max 


rd 


(2) 
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where 


Y » the  total  unit  weight  of  the  soil 


H 3 depth  from  the  ground  surface  to  the  point  in  question  in  feet 


g = acceleration  of  gravity 


the  peak  acceleration  at  the  ground  surface  generated  by  the 
earthquake  in  the  same  system  of  units  as  g 


a rigidity  factor 


The  constant  0.65  is  a factor  which  corrects  the  maximum  shear  stress 


to  an  equivalent  sinusoidal  shear  stress 


can  be  determined  for  any  depth  in 


Using  equation  2,  t 


factor  used  in  this  analysis  is  shown  by  the 


the  soil  profile.  The  r 


curve  marked  "analysis"  on  Figure  8 


Earthquake  parameters 


17.  In  order  to  conduct  this  analysis,  the  maximum  acceleration 


generated  by  the  design  earthquake  and  the  number  of  equivalent  cycles 


of  stress  are  required.  The  maximum  acceleration  was  specified  by  the 


NEC  as  0.12  and  0.2  g.  Review  of  the  geological  and  seismological 


studies  conducted  at  LaCross  predict  that  an  earthquake  of  Modified 
Mercalli  Intensity  VIII  in  the  epicentral  region  has  occurred  on  the 
Keweenaw  fault.1  For  analysis  purposes,  an  earthquake  with  an  intensity 
one  unit  greater  than  the  largest  recorded  intensity  was  assumed.  Thus, 
using  a Modified  Mercalli  Intensity  of  IX  and  using  the  intensity- 
magnitude  relationships  shown  on  Figure  9 (Reference  6),  a magnitude 
6.6  earthquake  is  postulated.  Figure  10  is  a plot  of  number  of  equivalent 
cycles  versus  magnitude  which  was  developed  by  Seed.  This  plot  was 
entered  and  10  equivalent  cycles,  which  are  essentially  an  upper  bound 


to  the  data  in  the  plot,  were  assumed  appropriate.  For  the  design 
earthquakes  (SSEa,  and  SSEb)  10  cycles  and  peak  ground  surface  accelera- 
tion values  of  0.12  and  0.2  g were  used,  respectively. 

Analysis 

18.  Because  of  the  high  N values  obtained  in  the  dense  sand  layer 
at  a depth  of  about  105  to  115  ft,  this  zone  is  predicted  to  remain  stable 
even  under  the  so-called  design  earthquakes.  There  are  no  other  data 
available  from  this  layer.  Tf  one  assumes  that  the  dynamic  soil  strength 
of  this  zone  was  the  same  as  that  Judged  appropriate  for  the  upper 
materials  by  WES  on  Figure  7,  then  liquefaction  might  be  predicted. 

WES  does  not  believe  this  will  happen. 

SSE,  0.12  g 

19*  Table  1 presents  the  information  needed  to  calculate  the 

average  shear  stress  (*av.e)  for  the  soil  profile  using  the  Seed-Idriss 
2 

method.  Also  listed  on  Table  1 is  the  effective  overburden  pressure 

(a  ) needed  to  enter  Figure  7 to  determine  the  available  soil  strength, 
o 

Both  the  values  of  dynamic  shear  strength  and  dynamic  shear  stress  are 
listed  as  a function  of  depth  on  Table  2.  The  factor  of  safety  against 
10  percent  double-amplitude  strain  has  been  defined  as  the  dynamic  shear 
strength  divided  by  the  average  dynamic  shear  stress.  This  factor  of 
safety  is  also  listed  on  Table  2.  It  should  be  noted  that  the  factor 
of  safety  below  a depth  of  35  ft  (depth  of  excavation)  varies  from 
0.99  to  1.15  and  is  below  1.1  to  a depth  of  100  ft.  Factors  of  safety 
less  than  1 were  predicted  in  the  soil  between  a depth  of  35  and  I4 5 ft. 

It  should  be  emphasized  that  the  state  of  knowledge  is  not  adequately 
refined  and  the  assumptions  required  to  carry  out  the  analysis,  given 
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the  limitations  in  the  existing  data  base,  are  such  that  it  is  not 


believed  Justified  to  call  this  site  safe  even  though  factors  of  safety 


marginally  greater  than  1 were  calculated.  However,  this  analysis  does 


not  prove  the  site  unsafe  under  this  acceleration  as  it  is  possible  that 


had  more  extensive  data  and  more  thorough  documentation  been  available 


the  Judgments  concerning  the  in-situ  density  and  cyclic  shear  strength 


would  have  been  different 


20.  The  information  needed  to  calculate  t 


for  the  SSE  with  0.2  g peak  acceleration  is  given  on  Table  3-  Also  shown 


on  Table  3 is  the  effective  overburden  pressure  which  was  used  to  enter 
Figure  7 to  determine  the  shear  strength  of  the  soil.  Table  ^ presents 


the  average  shear  stress  and  dynamic  shear  strength  as  a function  of 


depth  for  this  SSE.  Also  shown  on  Table  h are  factors  of  safety  (as 


previously  defined)  for  this  SSE.  As  can  be  seen,  the  factors  of  safety 
below  a depth  of  35  ft  vary  from  0.59  to  0.66.  Clearly,  this  indicates 


failure,  as  failure  is  defined  in  this  report.  A doubling  of  the  cyclic 


strength  over  that  shown  by  the  authors  would  be  required  to  produce  a 


This  level  is  often  considered  reasonable 


Empirical  Liquefaction  Analysis 


Empirical  data  in  the  form  of  stress  ratio  and  corrected  N 


values  for  sites  that  have  and  have  not  liquefied  during  past  earthquakes 
have  been  developed  and  plotted  on  Figure  11  (Reference  5)-  On  Figure  11 


the  Standard  Penetration  Test  N values  have  been  corrected  to  an  effective 


) 

overburden  pressure  of  1 tsf.  As  a second  means  to  evaluate  the  lique- 
faction potential  at  the  LaCross  Site,  the  average  blows  per  foot  as  shown 
on  Figure  1,  were  plotted  against  the  stress  ratio  and  compared  to  the 

data  shown  on  Figure  11.  The  stress  ratio  t divided  by  o where 
B ave  o 

is  equal  to  the  effective  overburden  pressure  is  also  tabulated 
on  Tables  1 and  3 for  the  various  depths  in  question. 

22.  Tiie  blows  per  foot  were  assumed  to  be  SPT  N values  and  were 
corrected  to  an  overburden  pressure  of  1 tsf  by  the  formula: 

- CN  x N (3) 

where 

N = Standard  Penetration  Test  penetration  resistance  value  measured 
in  the  field 

= Standard  Penetration  Test  N values  corrected  to  an  overburden 
pressure  of  1 tsf 

C.T  * Correction  factor 
N 

23.  Cjj  was  determined  from  Figure  12  (References  5 and  7). 

Figure  12  was  extrapolated  back  to  zero  using  data  in  Peck,  Hanson,  and 
7 

Thornburn.  Values  of  C„  and  N,  are  also  shown  on  Tables  1 and  3. 

N 1 

These  values  have  been  superimposed  on  Figure  11.  Most  of  the  values 
show  that  liquefaction  should  not  occur;  however,  liquefaction  is  predicted 
from  25  to  35  ft. 

2U . A similar  analysis  was  conducted  assuming  a SSE  with  a peak 
acceleration  of  0.2  g.  The  stress  ratios,  , N , and  , as  a 
function  of  depth  are  tabulated  on  Table  3.  These  values  are  shown  on 
Figure  13.  These  data  indicate  that  liquefaction  is  possible  if  a SSE 
producing  0.2  g at  the  ground  surface  occurs.  The  data  points  to  the  far 

lU 


i 


right  which  indicate  safe  conditions  are  for  the  10-ft-thick  sand  and 
gravel  layer  at  a depth  of  about  105  ft  where  the  pileB  end. 

25*  Based  on  the  Japanese  experience  in  the  16  June  196**  Niigata 
earthquake,  a rule  of  thumb  has  been  developed.  This  rule  states  that 
in  order  to  be  safe  against  liquefaction,  the  Standard  Penetration  Test 

Q 

N value  should  be  at  least  two  times  the  depth  in  meters.  A line 
has  been  drawn  on  Figure  2 indicating  what  the  N value  should  be  if  it 
were  greater  than  two  times  the  depth  in  meters.  Note  that  a large 
majority  of  the  blows  per  foot  fall  on  the  unsafe  side  of  this  line. 

This  is  particularly  important  in  this  case  since  the  peak  acceleration 
at  Niigata  was  approximately  0.l6  g. 

Summary  and  Conclusions 

26.  The  liquefaction  potential  of  the  LaCross  Site  was  evaluated 
for  two  earthquakes;  namely,  a SSE  with  a peak  acceleration  at  the  ground 
surface  of  0.12  g,  and  a SSE  with  a peak  acceleration  at  the  ground 
surface  of  0.2  g.  The  analysis  was  made  by  two  methods;  namely,  the 
Seed-Idriss  Simplified  Procedure  and  an  empirical  procedure.  Figure  lU 
is  a summary  plot  of  the  dynamic  shear  stress  as  a function  of  depth  for 
a peak  acceleration  of  0.12  g.  Superimposed  on  this  plot  is  a cyclic 
strength  of  the  material  assuming  10  equivalent  cycles  of  loading.  Note 
that  liquefaction  is  predicted  by  the  Seed-Tdriss  calculations  between 
a depth  of  32  and  U8  ft  and  liquefaction  is  predicted  by  the  empirical 
methods  between  a depth  of  2U  and  35  ft.  Japanese  experience 
at  Niigata,  Japan,  also  indicates  that  liquefaction  would  be  predicted 
below  a depth  of  15  ft.  As  indicated  in  Appendix  A,  the  piles  could 
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support  concentric  static  loads  even  after  a loss  of  lateral  support  down 
to  a depth  of  U8  ft.  However,  it  would  require  dynamic  structural  analysis 
heyond  the  scope  of  this  study  to  Judge  whether  they  would  have  sufficient 
bending  resistance  to  withstand  the  transient  eccentricity  of  the  static 
vertical  loading  and  the  transient  horizontal  loads  caused  by  the  seismic 
excitation.  In  view  of  the  fact  that  there  appears  to  be  no  reinforcing 
bars  in  the  top  one-third  of  the  pile 5 (as  called  for  in  some  seismic  design 
codes ^ ) , it  is  probable  that  the  available  bending  resistance  is  modest. 

27.  Figure  15  is  a summary  plot  of  the  dynamic  shear  stress  as  a 
function  of  depth  for  the  soil  profile  assuming  a peak  acceleration  of 
0.2  g.  Also  shown  on  this  plot  is  the  dynamic  shear  strength  of  the  soil 
assuming  10  equivalent  cycles  of  loading.  The  Seed-Tdriss  Simplified 
Procedure  predicts  liquefaction  below  a depth  of  25  ft.  The  empirical 
method  predicts  liquefaction  between  a depth  of  25  and  60  ft  and  between 
a depth  of  85  and  105  ft.  If  lateral  support  is  lost  in  the  depth  ranges 
predicted  by  either  method,  the  piles  would  be  in  danger  of  buckling 
failures  as  indicated  in  Appendix  A. 

28.  Based  on  the  Judgments  concerning  the  density  and  strength 
data  and  on  analysis  as  presented  herein,  the  soils  below  the  reactor 

at  the  LaCross  Site  are  predicted  to  strain  badly  if  a SSE  which  produces 
0.12  g at  the  surface  of  the  soil  occurs.  The  soils  beneath  the  reactor 
vessel  at  the  LaCross  Site  are  predicted  to  experience  excessive  strains 
and  liquefaction  if  the  SSE  with  a peak  acceleration  at  the  ground  surface 
of  0.2  g occurs.  Because  of  the  limitations  in  the  current  state  of 
knowledge  concerning  liquefaction  and  because  of  the  limited  data 
available  for  use  in  this  analysis,  WES  cannot  conclude  that  the  reactor 
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^KfSt 
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vessel  foundation  is  safe  if  the  0.12  g SSE  occurs  and  concludes  that 
the  reactor  vessel  foundation  is  unsafe  if  the  0.2  g SSE  occurs. 
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Table*  2 


RESULTS  OF  DYNAMIC  ANALYSIS  0.12  g SSE 


Depth 

ft 

Dynamic 

Shear 

Strength 

psf 

Average 

Shear 

Stress 

T 

ave 

Factor 

of 

Safety 

2.5 

50 

21 

2.38 

7*5 

102 

63 

1.62 

11.5 

1**5 

96 

1.51 

lU.O 

166 

116 

1.1*3 

17.5 

185 

li*6 

1.27 

22.5 

212 

188 

1.13 

27.5 

239 

229 

1.01* 

32.5 

293 

296 

0.99 

U2.5 

320 

321 

0.99 

1*7.5 

3>*7 

3l*2 

1.02 

52.5 

37U 

365 

1.02 

57.5 

1*00 

385 

1 ,0U 

62.5 

1*29 

1*11* 

1 .01* 

67.5 

1*60 

1*38 

1.05 

72.5 

1*90 

1*65 

1.05 

77.5 

520 

1*93 

1.05 

82.5 

551 

519 

1.06 

87.5 

581 

5l*l* 

1.07 

92.5 

611 

568 

1.08 

97.5 

61*2 

591 

1.09 

102.5 

672 

613 

1.10 

f 


I 

3 


« « 

a 3 , 

t ™ =*  * 

3 - 


£3*1 


g „ 

•H  O 
♦>  ♦» 
o o a 
4>  a o 

3 


5S  s 

h 0] 

v»  83  H I 
W CU 


SP  «d  • « vJ 

»«  « « >10 

V X U < 

> &5  4f  H 


nNHo<QS)«o>J  n<y(y|H>vo«oiniA44  n w 

r4  H H HHHHHHHHHHHHHH 


I g<aH 


33 


• ij 

53!  - 

4>  «>  « o < 
o jq  u o 
v a 4Jio  . 

4>  w 

S3 


3 S S 5 


cococOrHmvoi/\oinrHt—  -4  »Hao»/vcMOooin_4  mrHoa3t—voin-4 
-4  CM  rH  H O O Ov  O CO  03  h t-  t-  VO  VO  >5  (AIAlTklfVIA«\JJ^-IJ 

r-Ir^rHrHrHoddoddddooddodddoddoood 

(VJ  CO  tf)  (U  VO  H (M  VO  r«-  »/\  rH  cn  IA  IA  (\J  rH  O 0vt“V04(\iHfflV0t*lAlA 

n(V|(MM4VO(-ht-t-t-'i'0'OVOVOiOirtlrtlAtfV(Airi44444 
H H H H H H H H rH  H H H pH  H pH  H pH  rH  pH  pH  rH  rH  H H pH  pH  pH  H 

ooodoodooooooooooooooooooooo 

VOl/VO-4cnmCMO<’r>l/,VrHCOrHOC?vVOrH-4irvVO-4  r-n  co  t~-  o o t*  ao 
m O vO  CA  -4  rH  CO  -4  Ov  fO  t-  O -4  Ov  (VI  I-  N VO  0-400  CVJUNOCOt©  0-4 

OrHfHrH{\J0nr*l-4-4  »r\l/\VOVO\0*‘-t-a><OOvOvOvOOOrHrHCMCM 

ddddddddd^dddodddddddorHrHrH  rHrHrHrH 

mroropnrommfnrom^mmf^mmmrofnmmfnfnrommroro 

HHHHrHHHfHHHHHHHHHHHHHHHHHHHrHH 

ooodddododoooooooooooooooooo 

oooooooooooooooooooooooooooo 

CVICMCMCMCVlCMCMCMCgCMCMCMCVlCMCMCMCMCVICMCVjCMCMCMCMCMCMCMCM 

ooodoodooooooooooooooooooooo 


rHOOOOOOOOOOOOOOOO 


OOOOOOOOO 


StA\0(MHt(v>nio(noi4  rH  o o t-  vo  i/v  • m cu  h vo  m m cm  rH 

-4  mcOrHOVCOVOl.VrOCMO  O O Q O rH  CM  CO  -4  IT.  VO  CD  (VI  IA  « [■ 

H IA  CO  O rH  CVJ  -4  VO  <0  O CM  -4  in  r-  O CM  -4  VO  CO  O CM  -4  VO  OvrHrnin 


o o o o 


.H  CVI  CM  CM  <M  CM 


cnrocnro-4-4-4-4 


(A  IA  IA 


CM  lTv  VO  CO  O CM  I/V  Oi  CM  IA  Ov  CM  UV  rH  I-  -4  rH  t-  -4  H > 4 H H H •J’ 

04(nMC>9tcOh-|-VOIAIA4  inVOcOOrHrOmVOOOrHVOOrHO^  VO 
rHW>cOOCMVDO-4aiCM'00-4«OCM'5rHITvO\COt-rHVOOirvOv<nvO 


CM  CM  CM  mm-4_»-4  in  in  VO  VO  VO  t-  »—  CO 


Av  Ov  Ov  o O 


(a  o 4 cot-mrorHOii“-irifOrHf-o<'r*'o  th  w ia  o h ov  w o « vo  vo 

r-  H IA  h-  VO  4 CMOr-tfVrOrHOl<OOrHCMCOl/\'£t'-  C>>CMOvrn.4in  p 
<MCOCM_4VOavCMM>K-OfnM3cO  rH  IA  CO  rH  -4  t—  OrOVOOAAt—  OrOVO 

doHHHH('iNoifn(niA(A4  4 4 lAIAIAVOVOVO^h^COOOflO 


(A  CO  4 incOCQflOCO<OaOOOCOCpVO  rH  VO  rH  VO  rHVOrHVOVOrH  rHVOrHrH 

CSKK«52(;H3SSf;8S(5l:S8SS:SSS£S8 

ddHHHNMin4  4IA'0\0h>(0OodHHC'Ic'i(n4  4 IA'0 

•HrHrHrHrHrHrHrHrHrH 


IAIAIAOIAIAIAIAIAIAIAIAIAIAIAIAIAIAIAIA 
CM  N 


IA  IA  IA 


3iSS!Ut;2JS!KS£fifS£S;6;nS5»GK 

H-c.rrr»l..l»o3!!||;!r  r.h.JKjJD 


- **•- 


5.T37 


Table  k 


DYNAMIC  ANALYSIS  FOR  0.20  g SSE 


Depth 

ft 

Dynamic 

Shear 

Stress 

psf 

Average 

Shear 

Stress 

T 

ave 

Factor 

of 

Safety 

2.5 

50 

30 

1.39 

7.5 

102 

105 

0.97 

11.5 

1U5 

160 

0.91 

lU . o 

166 

191* 

0.86 

17.5 

185 

21*3 

0.76 

22.5 

212 

313 

0.68 

27.5 

239 

382 

0.63 

32.5 

266 

1*1*0 

0.60 

37.5 

293 

1*93 

0.59 

1*2.5 

320 

535 

0.60 

1*7-5 

3U7 

571 

0.6l 

52.5 

371* 

608 

0.6l 

57.5 

1*00 

61*1 

0.62 

62.5 

1*29 

690 

0.62 

67.5 

1*60 

729 

0.63 

72.5 

1*90 

776 

0.63 

77.5 

520 

821 

0.63 

82.5 

551 

861* 

0.61* 

87.5 

581 

905 

0.61* 

92.5 

6ll 

91*6 

0.65 

97.5 

61*2 

981* 

0.65 

102.5 

672 

1021 

0.6  6 

Idealized  Soil  Profile  and  Approximate  Average  Values  Assigned  to  Each  Layer  for  Analysis 
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Figure  2.  Blow  Counts  Presented  ns  a Function  of  Depth 
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Figure  7 • Superimposed  Dynamic  Shear  Strength  Versus  Normal  Effective 
Stress 
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Figure  8.  Rigidity  Factor  Versus  Depth  (after  Reference  2) 


Relation  Between  Earthquake  Magnitude  and  Intensity  in  Western  United  States 

(after  Reference  6) 
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Figure  10.  Equivalent  Numbers  of  Uniform  Stress  Cycles  Based 
on  Strongest  Components  of  Ground  Motion 
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Correlation  Betveen  Stress  Ratio  Causing 
Liquefaction  in  the  Field  and  Penetration 
Resistance  of  Sand 


Figure  11 
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Correlation  Between  Stress  Ratio  Causing 
Liquefaction  in  the  Field  and  Penetration 
Resistance  of  Sand 
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Figure  lU.  Results  of  SPT  Etapirical  Study  for  SSE  * 0.12  G 


APPENDIX  A 


CALCULATION  OF  BUCKLING  LOADS  FOR  PILES 

While  the  containment  vessel  foundation  mat  is  supported  on  piles 
on  a hearing  layer  which,  in  WES'  opinion,  will  not  liquefy,  these  piles 
require  lateral  support  above  this  layer  to  prevent  buckling  failure. 

If  liquefaction  occurs  in  some  depth  region  along  the  pile,  this  support 
will  be  dimenished  or  absent.  Further,  as  upward  seepage  from  the 
liquefied  zone  develops,  lateral  support  may  be  lost  all  the  way  from  the 
zone  to  the  base  of  the  mat  foundation.  The  Euler  buckling  load,  Pcr  , 
for  a typical  pile  has  been  calculated  as  follows: 

2.05ir2EI 
cr  “ L2 

where 

E = Young's  Modulus  of  the  pile 

I = Moment  of  inertia  of  the  cross  section  about  its  neutral  axis 

L = Unsupported  length 

The  following  assumptions  were  made: 

a.  The  pile  was  9 in.  in  diameter.  It  actually  had  an  8 in. 
diameter  at  the  tip  and  a 12  in.  diameter  at  the  butt. 

b.  The  pile  was  made  of  a linear  elastic  material  with  a 
Young's  Modulus  of  3,000,000  psi.  It  actually  was  made  of  3500  psi, 
28-day  strength,  cast  in-place  concrete  inside  a thin  steel  shell. 

c.  The  pile  was  fixed  at  the  base  of  the  mat  foundation. 

d.  The  pile  was  pinned  at  a depth,  L , below  the  base  of 

the  mat. 

For  these  assumptions,  the  relation  of  P to  L is  shown  on  Figure  Al. 

cr 


Al 


-w 


These  piles  have  been  rated  in  Reference  10  to  have  a 50-ton  static 
load  capacity.  Presumably,  the  piles  have  vertical  loads  considerably 
less  than  this  value.  At  the  rated  load,  the  analysis  indicates  that  the 
unsupported  length  at  which  buckling  would  take  place  is  approximately 
Uo  ft.  On  Figure  lU  of  the  main  text,  the  vertical  distance  from  the 
mat  to  the  bottom  of  the  shaded  zone  is  2U  ft  or  less.  The  piles  appear 
capable  of  supporting  their  vertical  working  loads  even  if  lateral  support 
is  lost  in  this  region.  However,  there  are  also  horizontal  dynamic  loads 
that  act  on  the  pile  butts  as  a result  of  earthquake  excitation  which 
tend  to  bend  the  piles  in  themselves  and  make  the  vertical  loads  eccentric 
which  would  in  turn  cause  further  bending.  Some  rough  calculations  suggest 
that  the  bending  capacity  of  the  piles  is  low  with  respect  to  the  moments 
which  might  occur.  A more  thorough  investigation  of  the  dynamic  bending 
problem  is  a complex  but  tractable  structural  dynamics  problem  beyond  the 
scope  of  this  study. 

For  the  0.2  g loading  there  is  no  point  in  performing  such  an  analysis 
as  Figure  15  of  the  main  text  indicates  a possible  unsupported  length  of 
over  80  ft.  If  loss  of  lateral  support  should  occur  over  this  length, 
as  shown  on  Figure  Al,  the  pile  would  buckle  under  its  static  load  alone. 
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In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 
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